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Abstract. The pair distribution function (PDF) is a key quantity for the analysis of correlation 
effects of a quantum system both in equilibrium and far from equilibrium. We derive an 
expression for the PDF in terms of the single-particle Green's functions — the solutions of the 
Keldysh/Kadanoff-Baym equations in the two-time plane — for a one- or two-component system. 
The result includes initial correlations and generalizes previous density matrix expressions from 
single-time quantum kinetic theory. Explicit expressions for the PDF are obtained in second 
Born approximation. 



1. Introduction 

In recent years the interest in correlated quantum many-body systems has increased steadily 
due to the observation of correlation effects such as liquid and crystal formation, superfluidity of 
Bose systems or bound states. Examples of systems include nuclear matter [TJ [2] , solid state 
systems |3] and dense astrophysical, laboratory dusty plasmas or the quark gluon plasma, for a 
recent overview see |H E] • A quantity sensitive to spatial correlations is the pair distribution 
function (PDF) h{r) — the probability to find an arbitrary pair of particles at a distance r. 

Of particular current interest is the short-time behavior of correlated systems following 
external perturbation such as excitation by intense radiation, e.g., In quantum systems 

such as atoms or solids this is often connected with rapid electron thermalization coupled to the 
dynamics of electronic correlations. Time and space-resolved measurement techniques detecting 
chemical reaction products, nuclear collision fragments or electrons and ions produced by laser 
ionization, e.g., (8], are becoming increasingly powerful. This brings the direct measurement of 
time-dependent pair correlations described by the nonequilibrium generalization of the PDF, 
h(r,t), within reach and increases the need for theoretical and computational tools that are able 
to predict h(r, t). 

From the theory side PDFs are routinely computed for nonideal classical systems using, e.g., 
Monte Carlo or molecular dynamics simulations, e.g., [4 . Extensions to quantum systems in 
equilibrium are available, e.g., in the frame of quantum Monte Carlo methods. However, there is 
still a high demand for accurate nonequilibrium simulations. Theoretical access to nonequilibrium 
PDFs is straightforwardly reached within density operator theory, e.g., [2J [TUJ [TTJ . An alternative 



approach to quantum many-body systems out of equilibrium is based on nonequilibrium Green's 
functions (NEGF). Apart from exact time-dependent numerical calculations, e.g., JT2J, which 
are limited to a few particles, numerical results for the two-time NEGF are expected to be the 
most accurate ones for strongly correlated quantum systems in nonequilibrium. This concerns, 
in particular, full two-time calculations, i.e. direct solutions of the Keldysh/Kadanoff-Baym 
equations (KBE) which, after the pioneering work of Danielewicz [T], has now become routine 
[21 [IHJ HH H51 ESI E3 EE]- They yield the single-particle Green's functions and all one-particle 
observables, including spectral function, as a function of time. A number of two-particle 
quantities can also be computed, taking advantage of the two-time structure of the equations, 
most importantly the interaction energy. However, in this approach the two-particle Green's 
function is eliminated by introduction of the selfenergy £ (the integration contour C and the 
arguments of the functions will be explained in Sec. [2|, 



This concept of the selfenergy has proven extremely successful for reducing the many-particle 
problem to an effective single-particle one (quasiparticle picture) which forms the basis of Green's 
function theory and the Feynman diagram technique. At the same time, the direct access to 
the pair distribution function — which contains important additional information beyond the one 
carried by the single-particle Green's function — is lost. Due to this reason, so far, no NEGF 
results for the full pair distribution function have been obtained which go beyond the computation 
of double occupancies, accessible from Eq. ([I]) 

In this paper we solve this problem. Starting from the Bethe-Salpeter equation for the 
two-particle Green's function, we demonstrate in Sec. [2] how, from a given time-dependent 
solution for the single-particle Green's function in a chosen approximation for the selfenergy 
S, the PDF in the same approximation can be reconstructed. In this paper, we choose, as an 
example, the second order Born selfenergy approximation E 2B which includes the interaction 
up to second order (E 2B oc V 2 ). Our results include the standard PDF h{r) in equilibrium, and 
its extension h{r\,T2) to inhomogeneous systems and to multicomponent systems where the 
PDF becomes a matrix /i ab (ri,r-2), a,b labeling the particle species. Further, we obtain results 
for an arbitrary nonequilibrium system where the PDF becomes time-dependent, /i ab — > h ah (t). 
Finally, the results are directly generalized to two-time PDF's, /i ab (ri,r-2) —> /i ab (ri, r-2, t). Our 
results selfconsistently include the influence of a correlated initial state, and they describe the 
decay of initial correlations which is particularly important at the initial stage of the relaxation. 
This problem was studied in detail for the single-time PDF within density operator theory 
[21 QUI EH 123) an d these results are contained in our theory as a special case. 

In Sec. [6] we discuss, as an illustration, the possible application of the PDF reconstruction 
algorithm to a two-component system of electrons and holes in a bilayer structure. This system 
has attracted substantial interest in recent years because it exhibits strong electron-electron and 
hole-hole correlations leading to spatial ordering and, at the same time, electron-hole bound 
states (excitons). There have been detailed investigations of the phase diagram |21[ [221 [231 I24| . 
of exciton and hole crystallization \2'2[ [231 126| . of collective excitations |27| |2"8] and of exciton 
Bose condensation and superfluidity, e.g., [251 122]- We conclude with a summary and discussion 
in Sec. 

2. Bethe-Salpeter equation 

In nonequilibrium Green's functions theory, the PDF of particle species "a" and "b" is obtained 
from the two-particle correlation function |5j taken at four equal times t\ = t2 = = t' 2 = t, 
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where and are fermionic or bosonic creation and annihiliation operators, and we introduced 
the short-hand notation 1 = Ti,s z i,t\. To determine the two-particle correlation function g 3 ^^ 
in nonequilibrium we start from the more general two-particle Green's function g 3 ^ on the 
Schwinger-Keldysh contour C which consists of a Hartree-Fock and a correlation part 

5 ab (12,l'2') = 5 | I b F (12,l'2') + < 7c a b rr (12,l / 2'), (3) 
S&(12,1'2') = 5 a (l,l')<7 b (2,2')±W(l,2')5 b (2,l'), (4) 

where we allow for four different time arguments, ti,t2,t'±,t 2 , all located on the contour C. The 
correlation part, <? ab rr , obeys the Bethe-Salpeter equation, e.g., 



<& b r (12, 1'2') = i J dld2dld2 <f(l, l)g h (2, 2)K ab (12, 12) 5 ab (12, l'2 ; ) , (5) 

which presents a formal closure of the second equation of the Martin-Schwinger hierarchy for 
the two-particle Green's function. Compared to Ref. |34j we also restored the exchange term 
where plus (minus) refers to bosons (fermions). In Eq. (|5| the function K ah is a general dynamic 
interaction (for details see Ref. [5]), which will be simplified in the following. 

1.) The first simplification is to introduce the screened ladder approximation, i<C ab (12, 12) — > 
V ab (12)S(l — 1)5(2 — 2). Suppressing for a moment the space integration and spin summation 
(they are implied by the repeated arguments 1 and 2 under the integral), we obtain from Eq. (pM) 



gt r (12, 1'2') = i J dhdt 2 g*(l, l)g h (2, 2> ab (12) 5 ab (12, l'2') . (6) 

We mention that general nonequilibrium initial correlations can be included by the proper 
definition of the contour C or accounted for via additional contributions to the single particle 
selfenergy, e.g., [TU EH 1321 133] • Equation ^ contains strong coupling and dynamical screening 
effects leading to a complicated integro-differential equation. 

2. ) Our next simplification is to neglect dynamic effects in the potential (dynamical screening, 
related to the GW approximation) which leads to the replacement l/ ab (12) — > V &h (fi2)5(ti — h), 
where fyi = | r i — r2 1 , 

5c ab r (12, 1'2') = i jf at g a (l, l)g b (2, 2> ab (f 12 ) 5 ab (12, l'2') , (7) 

and under the integral 1 1 = 1 2 = i. This equation corresponds to a static T-matrix approximation 
for the two-particle Green's function and the pair correlations, i.e. to a complete summation of 
the Born series. 

3. ) Since our goal is to reconstruct the pair distributions from a solution of the two-time 
KBE for the single-particle Green's functions in second order Born approximation, we limit 
ourselves to the first iteration of the integral equation (J7|), i.e. we replace, under the integral, the 
two-particle Green's function by the Hartree-Fock approximation: 

«f b (12, 1'2') =ijf di<7 a (l,lV(2, 2)F ab (r 12 )<£ b F (12, l'2') , with h = i 2 = t. (8) 
We can bring this equation to a more compact form by introducing the following new functions 
G ab (12,l'2') = 5 a (l,lV(2,2'), (9) 

S ab (12, 1'2') = y ab (r 12 ) {Gg b (12, l'2') ± 5 ah G$° '(12, 2'l') } , (10) 



where Gq 13 is the Hartree approximation for the two-particle Green's function (g^p without 
exchange terms) and Y,q° the first-order two-particle selfenergy. With these definitions Eq. ^ 
becomes 

gf ori (12, 1'2') = i / dt Gg b (12, 12)Sg b (12, l'2') , with t x = t 2 = t. (11) 



c 



We underline that all expressions, so far, are written on the Schwinger-Keldysh contour, i.e. all 
single-particle Green's functions are 3x3 matrices whereas the two-particle functions contain 3 4 
Keldysh matrix elements. To obtain the nonequilibrium pair correlation function we now have to 



extract from Eq. (11) functions depending on the real physical time. In particular, according 
to Eq. we need only the two-particle correlation function which we consider in the 

following. 

3. Reduction of the time structure of the two-particle correlation function y ab < 

In order to compute the pair distribution function, Eq. ([2]), the two-particle correlation function 
(7 ab,< at four equal times is needed. To this end, we first determine the equation of motion for 
the two-time two-particle correlation function g ab,< (i, t') and then specialize to the one-time 
two- particle correlation function g ab,< (i). 

3.1. Two-time two-particle correlation function g ^it^t') 

The problem of the large number of Keldysh matrix elements of g ab can be simplified drastically in 
the static Born approximation. As a first step in simplifying the time dependencies we specialize 
to functions with pairwise equal time arguments, t\ = i 2 = t, and t[ = t' 2 = t' . Then it is clear 



from Eq. (11 ) that the functions Gfj and T,^ 3 depend only on two times in the following way 



G^inrry^tt') = iftnty^fat-yj), (12) 

^(r^r'^W) = F ab (r 12 ){Gg b (r 1 r 2 ;r / 1 r' 2 ;tt / )±^bG ! g b (r 1 r 2 ;r / 2 r / 1 ;tt / )} , (13) 
where we also restored the coordinate arguments (spin variables are not written explicitly). Then 



Eq. ( 11 ) turns into 

g^r^r'^tt') = i J dt Gg b ( ri r 2 ; f x r 2 , tt)Sg b ( ri r 2 ; r'^; it') . (14) 

With this result, Eq. (|| for g ah can be rewritten as 

5 ab ( ri r 2 ; rir'a; tt') = <? H b F (rir 2 ; r'^; tt') + <7 c ab rr (rir 2 ; r^; tt') , (15) 



ab 

the form 



where g^ VI is completed by the Hartree-Fock contribution which, in the new notation, acquires 



< F ( ri r 2 ; r' ir ' 2 ; tt') = G ab (nr 2 ; r[r' 2 ; tt') ± 5 ab Gg b (nr 2 ; r' 2 r[; tt') . (16) 



According to Eq. (J2J), g ab (ii') reads in terms of creation and annihilation operators (suppressing 
spatial arguments), 



9 



*\tt')=i\*i(t')*l(t')* h (m,(t)). (17) 



It is important to notice that this structure is different from a generalized density-density 
correlation function 

cl b n (tt') = ?(n h {t')n{(t)) = i 2 (*l(t')* h (t')*lm.(t)) + <f b (if) , (18) 



since the order of the operators, in general, cannot be interchanged. 



Now, to compute the pair distribution function from Eq. (15), we have to extract g ab ' < from 
the Keldysh matrix and, in particular, the "<" component from the contour integral in Eq. (14) . 
The solution of this problem is well known for one-particle functions, and for two-particle functions 
it has been solved in Ref. [31]. However, the latter results are not needed here. Indeed, although 
we are dealing with two-particle functions Gq 3 and £ ab , in the present approximation, they have 
the same time-dependence (they depend just on t, t') as occurs in the case of single-particle 
functions in the integral term of the Keldysh Kadanoff-Baym equations. As a consequence, both 
these functions are 3x3 Keldysh matrices and, therefore, also the two-particle Green's function 



<7 ab , Eq. ([14]), is a 3 x 3 Keldysh matrix. 

The Keldysh components are classified by location of the two time-arguments on the contour 
C. They include correlation functions, ">, <" with two real arguments, "[,]" having one real 
and one imaginary time argument and the Matsubara component, a = M, where both arguments 
lie on the imaginary track. Below we will need only the "<" component which determines the 
time-dependent nonequilibrium PDF and the Matsubara component which yields the PDF in 
thermodynamic equilibrium. We start the analysis with the former and consider the latter in 
SecHJ 

Each of the nine components G ab ' a of the matrix Gq 3 is of the product form 

Gf'^nrajr^ttO = g^t; r'^g^^t; r' 2 t') , (19) 



which allows for an easy matrix multiplication of Gq 13 and £ ab under the integral in Eq. ( 14 ) 



Since each component a has the same dependence on the coordinates as in Eq. (14), we will 
suppress the space arguments and integrations in the remainder of this section. 

The result for the "<" component is well known (Langreth rules) and consists of an initial 
correlation and a collision term 

g^(t, t>) = g$<(t, t>) + g%<(t, t>) , (20) 

which arise, respectively, from the imaginary and real part of the contour C, 

g$<(t,t') = -i f"drG^(t,r)^ h ' l (r,t'), (21) 
Jo 

ottW) = jTdf{G^ < (t 1 ^E^(t,tO + Gf^(t 1 ^Ef , ' < (t > 0}. (22) 

Here GQ b,< denotes the product of two single-particle correlation functions <? a,< ,g b,< , both 
having two real time arguments, whereas the components ] and [ describe products of functions 
depending on one real and one complex (r) time argument describing the propagation of initial 
correlations, for a detailed discussion, see, e.g., Refs. [IUH7]. Further, we introduced retarded 
and advanced functions defined by 
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(t,t') = ±e[±(t-t , )]{a > (t,t')-a < (t,t')} . (23) 



3.2. One-time two-particle correlation function g ^(t) 

We now perform the final step on the way to the nonequilibrium PDF — taking the time-diagonal 
limit, t = t', in the above equations for which has the general structure 

5 ab '< W = + 5i a c < (*) + gSf® ■ (24) 



Expressing the time-diagonal part of the single-particle correlation function by the single-particle 
density matrix p a (i) = ±i<7 a (t, t), where the upper (lower) sign refers to bosons (fermions), we 
obtain for the two-particle correlation function on the time-diagonal 



G^<(t,t) = -p a (t)p h (t). 



(25) 



With this relation and Eq. (16) we first obtain the nonequilibrium Hartree-Fock pair distribution 
function 



^ < (rir 2 ;r / ir' 2 ;t) = -p^r 1 r' 1 ,t)p b (r 2 r' 2 ,t) T 5 ah p a (r 1 r' 2 ,t)p b (r 2 r' 1 ,t) . 



„a/„ _/ j.\ „b/„ „/ 



(26) 



The second contribution to Eq. (24) due to the initial correlations follows from Eq. (21 ) and is 
given by 



g^i^rW^t) = -i Jdrj d^d 3 ^^)^ (r x t; r x r)^ (r 2 t; r 2 r) 



(27) 



g^(v !T, ri, % b ' r (f 2 r, r' 2 ; t) ± W r (riT, r' 2 , t) 5 b ' r (f 2 r, r' i; t) I . 



Finally, the scattering contribution (22) contains two real-time integrals extending to t and t' , 



respectively. Since now both times are equal, a partial cancellation is possible. Indeed, taking 



advantage of relation (23) allows to cancel all contributions with products of four ">" or four 



: <" functions, and we obtain 



g^in^r'^t) =i |dtYdWr 2 F ab (f 12 )|g a ' > (r 1 t;f 1 l) 5 b ' > (r 2 t;r 2 t) (28) 

o 

x [g^irxt, ri, t)g h '<(f 2 t, r 2 ; t) ± 5 ab5 a '<( ri t, r 2 , t)g b '<(f 2 t, r' i; t)] 

This expression is readily understood: the first term describes the correlation build-up due to 
scattering of a particle pair a, b out of state | r i ) |f 2 ) into state |ri) |r 2 ), whereas the second 
contribution which is obtained by interchanging functions with the indices "<" and ">" describes 
the opposite process. The integral has the familiar non-Markovian structure ("memory") 
indicating that scattering processes from all times prior to the current one do contribute although 
their weight decreases since the single-particle Green's function decay with increasing difference 



of their time arguments. Expression (28) has exactly the same structure as known from density 



operator theory |llj . The binary density operator in Born approximation is recovered if the 
two-time correlation functions are eliminated using the generalized Kadanoff-Baym ansatz |37| . 
see also Ref. [35]. The present result is more general because this reconstruction ansatz, which is 
only valid approximately and becomes increasingly inaccurate with increasing correlation effects, 
is avoided. 



4. Nonequilibrium pair distribution function 

The nonequilibrium pair distribution function /j ab (t) follows immediately from the two-particle 
correlation function according to Eq. ^ and the results (24), (26), (27) and (28). Explicit results 



depend on the choice of a basis. We start from the coordinate representation. 



4-1- Coordinate representation 

To obtain the standard pair distribution function which is defined in configuration space we set 
ri = r'j and r 2 = r' 2 in all expressions 

^ b ( ri ,r 2 ,t) =p*(r 1 r 1 ,t)p h (r 2 r 2 ,t)±5 ah p a (r 1 r 2 ,t)p h (r 2 r 1 ,t) 

ab,</ ,\ ab,</ ,\ \ ' 

~9ic (rir 2 ;rir 2 ;t) - g co { (rir 2 ; rir 2 ; t) . 

In a spatially inhomogeneous system it is advantageous to introduce center of mass and relative 
coordinates, R = (ri + r 2 )/2 and r = ri — r 2 , leading to the replacements ri j2 = R ± r/2, on 
the right hand side. Then /i ab (R, r) is understood as the local probability to find a particle pair 
at distance r around space point R. 

If spatial inhomogeneities are not relevant or not of interest, a space integration leads to the 
space averaged nonequilibrium pair distribution 

h ah (v, t) = J d d R /i ab (R, r, t) , (30) 

where d is the dimensionality of the system. Note that /i ab is normalized to the total number of 
particles according to 

"d d h ah (r,t) = N & N h . (31) 



Finally, in some cases just the dependence on the magnitude of the pair separation is of interest 
which is obtained by an angle integration. The result is called radial distribution function, 
fo ab (r, t), and follows from Eq. (31) according to (we consider a two-dimensional system and 
change to polar coordinates) 



r2ir 

h ah (r,t)= dc/>r/i ab (r,<M), (32) 
Jo 

drh ah (r,t)=N a N h , (33) 



with the normalization 



and similarly in 3D. 



4-2. Arbitrary basis 

Let us now consider the case of an arbitrary orthonormal stationary basis {^(r)} which is of 
relevance, in particular, for spatially inhomogenenous systems. Then the single-particle Green's 
functions become matrices according to 

5^(1,1') =J2Mri)9ij a (h,t' 1 )ti(r' 1 ) , (34) 

ij 

which holds for any Keldysh component "a". Analogously, the two-particle single-time correlation 
function is represented by a four-dimensional matrix 

g^W^t) =Y,^)^)9^mUr[)^(r' 2 ) • (35) 

ijkl 



To compute the nonequilibrium PDF (24) we have to expand the Hartree-Fock, the initial 
correlation and the scattering part in this basis. For the Hartree-Fock contribution we obtain, in 
analogy to ( 35 ) , 

g^im^t) = E^( r i)^( r 2)^F^(*)^W)0r(r2), (36) 

ijkl 

&%U® = sJil'W^Wi^^W- (37) 



Similarly, we obtain from Eq. (27) for the initial correlation contribution 

ab,</ / / j\ 

Sic ( r i r 2;rir 2 ;i) 



ab,< /,\ 
9lC,ijkl\ l ) 



and from Eq. (28) 



ab,</ / / ,\ 

5coi (rir 2 ;r ir2 ;t) 



ijkl 
mnrs 

{£i(T, t)g h J(r, t) ± 5 ah £/(r, t) g% (r, t) } , 



£ ^(n)^(r 2 ) 5c %. w (i)^(r'i)0r(r' 2 ) , 



(38) 
(39) 



(40) 



q mnrs 

#r (f, *) ^ ft ± ^ab <;< (t , t)g b k >< (t, t) 
-(>**< and 5i i2 <H> #2,l) j • 



(41) 



Basis expansions for inhomogeneous quantum many-body systems have been successfully applied 
to electrons in quantum dots ("artificial atoms"), e.g., [35] . and small atoms and molecules, 
e.g., j!6] 136] . For these systems, the Keldysh/Kadanoff-Baym equations are solved for the 
matrix function gf-{t,t'). Using these results and formulas (36), (38) and (40), the two-particle 



correlation function, g> ab,< in configuration space, Eq. (35), can be reconstructed. Besides, also 
the matrix elements g^i (t) themselves are of interest, as they carry extensive information on 

the many-body system. For example, the matrix components 9ijij^(t) describe the correlation of 
two particles "a,b" occupying the orbitals i and j, respectively, at a given moment t. 

5. Equilibrium pair distribution function 

To compute the pair distribution in thermodynamic equilibrium from the Matsubara Green's 



function we return to Eq. (14) for the Keldysh matrix and extract the Matsubara component. 



According to Eq. (19) it is given by a product of one-particle Matsubara Green's functions 

<< M ( ri r 2 ; rW 2 ; r) = ^(nr', r) 5 b > M (r 2 r 2 , r) , 



(42) 



where r = t — t' . Similarly as in nonequilibrium, the two-particle Matsubara Green's function 
consists of a Hartree-Fock and correlation part, 



,ab,M/ 



9 (V = 9hf (V + Scoi ( T ) ( 43 ) 
where the latter is obtained from Eq. ( |19[ ) using the Langreth rules, in analogy to Eq. ( |21[ ), 



ab,M/ 



ab.M/ \ 

9coi ( T ) 



-l / drG (r-r)S (r) . 
Jo 



(44) 



Restoring the coordinate arguments and using the definitions of Gg b and £§ b , Eqs. g and §L0$, 
we obtain 

^ M (r ir2 ;r;r' 2 ) = -p^v'^v^) T «5 ab p a (r 1 r' 2 )p b (r 2 r / 1 ) , (45) 
^(rirasrir^T) = -\ J^df J d?nd?f 2 V^{f l2 ) x 

{^(nn.r-rj^.r-r) x (46) 
[^• M (r ir i, r) 5 b ' M (f 2 r 2 , r) ± W M (r ir ' 2 , f) 5 b ' M (f 2 ri, f )] } . 
From this the equilibrium PDF is obtained by introducing center of mass and relative coordinates 



as was done in the nonequilibrium situation, cf. Sec. 4.1 



6. Numerical example: pair correlations in an electron-hole bilayer system 

To illustrate the results obtained so far we consider an example of a two-component system 
where strong correlations play a prominent role. The system of interest consists of two layers 
of zero thickness which contain an equal finite number N of negative (electrons) and positive 
(holes) charged particles with a finite layer spacing d. In each plane, the particles are confined by 
a harmonic potential of frequency Vt. Such electron-hole bilayers have been actively studied in 
recent years, e.g., [2D E21 E3] ; because they allow to study strongly correlated excitons which 
may exhibit Bose condensation and superfluidity, as well as liquid-like and crystal-like behavior. 
For more details, we refer to Refs. 



6.1. Model 

The Hamiltonian of the quasi-two-dimensional electron-hole bilayer is given by 

H = H c + H h + H ch (47) 
ft - E ^ (-A,. + ■!.)+ A £ (48) 

1=1 Kj=2 y {Ti,e — Tj,e) 

N c N h 1 

#eh = -AEE 7 ==_=, (50) 

i=i j=i y' (r i)e - r j)h ) z + d* z 
where we introduced dimensionless variables by rescaling length and energy by the units, 

Eo = m - (51) 



The term m*^ h denotes the effective mass of electrons/holes and d* is the effective distance 
between the layers. Further, we introduced the coupling parameter A measuring the strength of 
the Coulomb interaction relative to the confinement energy 



Here a-Q is the Bohr radius of a Hydrogen-like bound state — an exciton (we use the electron mass 
instead of the reduced mass). 

The coupling parameter measures the strength of the Coulomb interaction among identical 
particles in each layer as well as the correlation between electron and holes. Generally, one may 
expect that for A < 1, i.e. for very strong confinement, there is a strong wave function overlap, 
and the system approaches ideal gas like behavior. In the opposite case, A>1, the Coulomb 
interaction dominates, and particles will tend to become localized. Finally, variation of the layer 
separation gives an additional control of the many-particle state: for d — > oo both layers will 
be decoupled, containing independent electron and hole populations, whereas for decreasing d 
Coulomb attraction plays an increasing role. This gives rise to formation of indirect excitons 
which behave (approximately) as bosons and exhibit dipole interaction, e.g., [26, 29 . 



6.2. Equilibrium PDF of the electron-hole bilayer — Comparison to Path Integral Monte-Carlo 
results 

Preliminary results for the equilibrium PDF of the electron- hole bilayer according to Sec. [5] were 
obtained recently [40 , however they still require further numerical tests. Therefore, to illustrate 
the physical content of the PDF, in this section, we show some results obtained from path integral 
Monte-Carlo (PIMC) calculations by Boning et al. . Their approach to the description of the 
bilayer system is different to the one presented in Sec. [6j as they assume beforehand, that the 
layer separation d and the interaction strength A in Eq. (47) between the electrons and holes, 
respectively, induces the formation of indirect excitons, i.e. quasipar tides, comprised of strongly 
bound but spatially separated electron-hole pairs. The excitons exhibit an interaction, which for 
large distances is of dipole type and for small distances approaches a soft Coulomb potential, 
for details see Ref. |29j . With these assumptions Boning et al. computed the equilibrium PDF 
/i XcX (r) of one exciton X relative to a fixed exciton X c in the center for different values of the 
exciton density n in a ZnS^Sei-^/ZnSe quantum well with doping factor x = 0.3. The density 
is measured in units of a^ 2 = 1.06 • 10 17 m -2 , where ag = ^ 2 e/(e 2 m*) is the electronic Bohr 
radius with the material constants e = 8.7 and m* e = 0.15 mo- The temperature is chosen to be 
kv,T = 0.001 Ha*, where the energy uni t is defined as Ha* = e 2 /(ea^) = 53.93 meV, resulting 



in a temperature of T = 0.63 K. In Fig. 6.2 /i XcX (r) is shown for different densities. One can 
see that for low densities, the excitons in the system are melted, showing no localization or 
radial ordering. If the density is increased, an exciton crystal starts to form, due to pressure 
crystallization. For further increased density this exciton crystal melts again, which is an effect 
of the Coulomb-like character of the exciton-exciton interaction at short distances, the excitons 
undergo quantum melting. The advantage of the PDF is clear from this figure: different phases 
of the system can be clearly distinguished which is not possible on the basis of single-particle 
quantities such as the density. 

PIMC is very efficient for computing the thermodynamic properties and also spectral properties 
[41 j of bosons, such as excitons. However, at high densities, excitons break up and form an 
electron-hole plasma (Mott effect) consisting of fermions. PIMC simulations of fermions at 
low temperatures are still hampered by the notorious sign problem |42j . In contrast, with 
nonequilibrium Green's functions, this regime is easily accessible whereas limitations arise 
with increasing coupling strength. Therefore, PIMC and NEGF have complementary areas of 
applicability. Moreover, NEGF should allow one to access the time-evolution of the PDF, as 
demonstrated above. 



7. Discussion 

In this paper we presented an approach to the pair distribution function of a quantum many- 
body system in the frame of the nonequilibrium real-time Green's functions. This problem is 
complicated due to the fact that the standard approach used in NEGFs uses a formal decoupling 




Figure 1. Radial Exciton-Exciton PDF h ° (r) in a ZnSo.3Seo.7/ZnSe quantum well at 
T = 0.63 K with layer separation d = 40.83 nm. Densities in units of 1.06 • 10 17 m~ 2 : (a) 
0.84 • 10" 3 , (b) 1.3 • 10~ 3 , (c) 1.7 • 10~ 3 , (d) 3.2 • 10~ 3 , (e) 3.6 • 10" 3 , (f) 4.0 • 10~ 3 . The radial 
length is measured in units of the so called Brueckner parameter r s = 3.26 • 10 s x a, where a is 
the mean interparticle distance. The illustration is taken from Ref. |29j . 

of the Martin-Schwinger hierarchy on the level of the first equation: the two-particle Green's 
function is eliminated by introduction of the single-particle selfenergy, cf. Eq. ([!]). With this 
elimination also direct access to the pair correlations is lost. Thus one has to reconstruct the 
pair correlations form the single-particle Green's function within a chosen approximation for the 
selfenergy. 

To solve this reconstruction problem we started the analysis from the equation of motion of 
the two-particle Green's function — the Bethe-Salpeter equation and simplified it systematically. 
We concentrated on the case of the static second Born approximation because for it a large 
number of numerical solutions of the KBE exist, for which it would be desirable to evaluate the 
PDF. It was shown that in second Born approximation a closed expression for the nonequilibrium 
PDF can be derived which can be straightforwardly evaluated. The result involves combinations 
of four single-particle Green's functions and is computationally expensive. Numerical results will 
be presented in a forthcoming paper. 
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